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We discuss possibility of describing solar, atmospheric and LSND resulta 
with four neutrinos forming two quasi-degenerate pairs. The simplest versions 
of this 2 + 2 scheme with either Ve or ly^ mixing exclusively with sterile neu- 
trino is disfavored by the SNO and atmospheric neutrino results respectively. 
A generalized scheme with sterile state participating in both the solar and 
atmospheric oscillations is stili allowed. We show that the complex pattern of 
mixing needed for this purpose follows from a simple Le+L^ — Lr — Lg symme- 
try. Specific form of Le + L^ — — Lg symmetric mass matrix is determined 
from experimental results. Two theoretical schemes which lead to this form 
and a proper breaking of + — L^- — Lg symmetry are discussed. 



I. INTRODUCTION 



Both positive and negative results on neutrino oscillation searches have provided very 
important clues on possible patterns of neutrino masses and mixing [|l|]. These results include 
experiments detecting solar and atmospheric neutrinos as well as laboratory experiments 
such as LSND, CHOOZ etc. The observed solar and atmospheric neutrino deficits provide 
concrete ground to believe in neutrino oscillations. Along with the LSND results, they give 
important and by now well-known |]T| Information on neutrino masses. 

The latest input in the analysis of neutrino spectrum is results of solar neutrino deficit 
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seen at SNO p[. This experiment finds lower neutrino flux in their charged current events 
compared to the flux inferred from the elastic scattering at SuperKamioka which receives 
contributions from the charged as well as the neutral current processes. The difference in 
these two fluxes is consistent with predominantly converting to active flavours only. The 
global analysis of the solar neutrino data ||^ shows that complete conversion of the solar Uf, 
to sterile neutrino is a disfavored possibility allowed at Sa level that too mainly in case of 
the vacuum oscillations only. 

The implications of the above results for the neutrino masses are as follows. 

• The oscillation interpretation of data requires IQ three different mass scales {/S.lsnd ~ 
eV^, Ayi ~ (5 — 8) ■ 10^^ eV^ and A5 < 10^'^ eV^) and four light neutrinos to account 
for these (mass)^ differences. 

• The schemes with hierarchical masses is highly disfavored from the point of view of 
explaining ali neutrino anomalies In this scheme, the allowed LSND probabili ty 
is found to be smaller than the observed one when negative results of neutrino 
oscillation searches at Bugey and CDHS are taken into consideration. 

• AH the experimental results before SNO could be understood ^ in terms of a simple 
picture in which four neutrinos group themselves into two pairs with a gap of order 



y/AisND- There are two versions: 

(a) in which one of the pairs consists of u^, Vs and accounts for the solar neutrino 
anomaly. The other is responsible for — Vj. oscillations of atmospheric neutrinos. 

(b) corresponds to converse possibility with — Vg accounting for the atmospheric 
neutrino deficit, solar neutrinos converting themselves completely to active compo- 
nents. Both these possibilities have been termed as 2 + 2 schemes of neutrino 
masses. 

The above possibilities can be incorporated naturally into theoretical schemes for neu- 
trino masses and there have been number of models [p|-[l3[] realizing these possibilities. The 
possibility (b) was already disfavored by the observed absence |T^ of the matter effects in 



atmospheric neutrino data. Now the SNO results ^ strongly disfavors (a). Thus it becomes 
a challenging task to understand experimental results within the four neutrino scenario and 
build necessary framework to account for this. We discuss simple four neutrino schemes in 
this paper which can simultaneously explain the solar, atmospheric and LSND results. 

While the possibilities (a) and (b) above are the simplest realization of the 2+2 schemes 
they are not exhaustive. In general, a sterile state can simultaneously but partially infiuence 
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both the solar and atmospheric neutrino oscillations. The most general possibility in this 
context was worked out in It was shown that the LSND results can be accounted for in 
2+2 scheme without conflicting with the negative searches of neutrino oscillations if v^, and 
reside mainly in different mass pairs. The unitarity of mixing matrix then dictates the 
foUowing mixing pattern among four neutrinos f^: 

= co^Qs v\ + sin 6^5' z/2 + C(e) , 
= cos Qa vj, + sin Qa Vi + , 
Vr = sin a(— sin 6 a 1^3 + cos 9a t-'a) + cos a(— sin 63 vi + cos 63 1^2) + , 

= cosq;(— sin 9a 1^3 + cos6'a ^'4) — sinQ;(— sin 63 + cosé'^ U2) + C(e) . (1) 



The v\^2 and ^3^4 represent two quasi-Dirac pairs of neutrinos with definite masses rriy^ {% = 
1, 4). 6s, Oa respectively denote mixing angles relevant for the solar and atmospheric neutrino 
oscillations. The angle a determines [| the amount of sterile component in these two fluxes. 
The splittings A12, A34 (Ajj = m^^, —mf^J respectively account for the solar and atmospheric 
oscillations. As is seen, the Ue and z/^ reside in different mass pairs apart from small 0{e) 
corrections. The parameter e thus decides the amplitude for the LSND oscillations. 

Two examples of pure active sterile mixing discussed above correspond to special cases 
of eq.(|l|) with a = tt/2 (case (a)) and a = (case (b)). While these extreme cases are ruled 
out, intermediate possibility with non- zero a is stili allowed. The phenomenology of this 
case was studied in [p|,p!5| and was updated after SNO results in [jl6[. Basically, the same 
parameter a determines amount of sterile component in the solar as well as atmospheric 
neutrino flux and thus can be constrained by both experiments. It was found in that 



the solar as well as atmospheric data can be fitted with a non-zero a in two possible ways. 
Either a is large in which case the sterile state mainly appears in the solar neutrino flux. 
The best fit in this case is obtained for 

sin^ a ^ 0.8 . (2) 

The possibility (a) above is allowed in this case only at 99%CL. The other case corresponds 
to small a and sterile component mainly in the atmospheric neutrino flux. The best fit value 
of a is given in this case by 

sin^a^O.l. (3) 



defined here differs from the one in by 7r/2. 
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The possibility (b) above can be allowed around 95%CL in this case. It follows that 
complete conversion of Ve to in the solar fiux or io Ug in the atmospheric fiux is a 
disfavored possibility although best fit vaine of a is not very far from these two limiting 
cases. General mixing pattern in eq.([^) is a more favored possibility. A priori, this mixing 
pattern looks complex but as we discuss below, it follows if 4X4 neutrino mass matrix 
respects a Lg + — — symmetry. We write down a specifìc ansatz based on this 
symmetry in the next section where we also show that + — Lr — Lg symmetric mass 
matrix can be determined from the experimental results. Next two sections are devoted to 
realization of the mass matrix based on Lg + L^ — L,- — symmetry. Last section summarizes 
the results obtained. 



II. Le + L^-Lr- Ls SYMMETRY AND NEUTRINO MIXING 



Let US parameterize the most general 4X4 matrix in the basis (z/g, i/^, z/^, z/^) as follows: 

(4) 



in^ A' 



Here each sub-block is a 2X2 matrix. The existence of pseudo-Dirac pairs implies some 
partially broken U{1) symmetry. Such symmetry follows if 2X2 matrices A, A' are sub- 
dominant compared to m in eq.(^. Let us choose A = A' = as a first approximation. 
Then for arbitrary m, A^j, is invariant under + — Lg — Lj. symmetry. It is possible to 
exactly diagonalize eq.(^ in this case. This is done in two steps. First, m is diagonalized 
by the foUowing 2X2 rotations R: 



Rfs m = Diag.{mi, 1712) 



(5) 



where a,f3 denote the angles of rotations and mi^2 are the eigenvalues of m. Given eq.(^), 
the foUowing 4X4 matrix diagonalizes ÀAi, when A, A' are zero: 



R'RAi^R^R'^ = Diag.{mi, —mi, m2, — ^2) 



(6) 



where, 



R 



( \ 
Rp ^ 

R^ 



(7) 
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and 



R' 



\ 

Eqs.(|^,|D together imply the following mixing among neutrinos: 

Ua = UaiVi = {R' R)li Vi . 
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~72 ) 



(8) 



(9) 



Explicitly, 



Vp 



v. 



-^(Z/I - U2) + - Z/4) , 

Sa / \ , Coi , N 

~V2 ~ ~ 



(10) 



The above pattern coincides with eq.(|I|) if S/3 <^ 1 and 9a = Os = ^- It is easily seen 
that /3 governs the amphtude of LSND oscillations. Eq . (p!0D imphes 

W 2eLSND ^ 4:{UeAU^4 + Ue^U^s f = sin^ 2f3 

when 77?,;^^ 2 -C m^^^ ~ eV. The observations at LSND then imphes that sin^ 2/3 ~ 3 ■ 10"'^. 

= 9s=Ti/^ is a prediction of the model which arise as a consequence of the + — 
Lj. — Lg symmetry. The perturbation A, A' which would cause the sphtting of the neutrino 
masses would also change this prediction somewhat but one will get two large mixing angles 
needed on phenomenological grounds. 

Eq.([TIl|) leads to the following survival probabihty in disappearance experiment such as 
CHOOZ [0] 

A,4t 



2 



4^--sin^2/3 



(11) 



This probability is correlated with the LSND angle. Its oscillatory part is suppressed in 
view of the LSND results on S/3. The average term is smaller than the present limit set by 



CHOOZ [1/] experiment but it can be signifìcant. 
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It is clear that the Le + L^ — Lr — Lg symmetric mass matrix leads to phenomenologically 
desirable pattern naturally. Its breaking is needed to generate the mass sphttings. Details 
depend upon the form and strength of the perturbation A, A' which is yet unspecified. 
But the following observation is relevant. One can calculate the (mass)^ differences using 
perturbation theory assuming that A, A' have general structure with ali elements having 
equal strength. Using the unperturbed eigenfunctions following from U in eq.(^, we get 

/^s = {ml^-ml;)^m^O{5) . (12) 

It follows that 

^ « ^O(l) , (13) 

A^ 1712 

Thus, hierarchy in the solar and atmospheric scales is linked in this case to intra-splitting 
between two pairs and one generically needs mi -C m2- We shall present an example where 
eq.(|T3D is realized naturally with 0{1) parameter being exactly 1. A^j^ is characterized 
by four parameters mi and m2- Ali these parameters are approximately determined 
phenomenologically: m2 and (3 from LSND results, a from general fit ||l^ to solar and 



atmospheric data and mi from eq . (p!3D . Thus M.^ gets phenomenologically determined in 
the symmetric limit. 

Let US parameterize m as follows: 



Oi a2 



m 



(14) 



\Ai A2 j 

The elements a, and Ai {i = 1,2) are determined in terms of mixing and masses as 
follows: 

A2 = m2Cf3Ca + rriiSpSa ~ m2Ca , 
Al = m2SaC(3 - miS/sCa ~ m2Sa , 
a2 = m2CaSi3 - miSaCp , 

ai = m2SpSa + rriiCfsCa ■ (15) 

The parameters on the RHS are directly determined from experiments. Since both sp 
and mi are required to be small, the above equations imply that 



and also 

which itself is small. Successful model should realize the above hierarchy and we present 
two specific examples. 



III. RADIATIVE SCHEME 

We showed that Le + L^j, — Lr — Lg symmetry leads to phenomenologically consistent 
2 + 2 model and identified a specific structure for the four neutrino mass matrix in this limit. 
We now discuss a radiative scheme which realizes this structure and also provides necessary 
breaking of the + — — Lg symmetry. Our starting point is the observation that a 
part of Le + Lfj^ — Lt- — symmetry namely, Lg + — Lj. already arise as an accidental and 
approximate symmetry in the Zee model |]l8l of neutrino masses. Thus it is naturai to start 



with this model. We add a singlet neutrino Vs to it and impose an associated singlet lepton 
number symmetry U{l)s which is carried by Vg-, rest of the fermions remaining unchanged 
under it. Due to inherent L^ + L^ — symmetry of the Zee model, the radiatively generated 
4X4 matrix automatically displays approximate Le + — L,- — Lg symmetry although we 
start with only a U{l)s symmetry. The Higgs fìelds of the model are two doublets 0i,2 and 



a charged singlet hi as in Zee model |18| and an additional charged singlet h2 carrying the 



same U{l)s charge as Vg. The leptonic Yukawa couplings in the model are then given by 

- Cy = fijllL ^jL K + PiUgCiR + QijìiL ejR(pi + H.C. . (18) 

Note that the U{l)s symmetry forbids the Dirac coupling between ii and Ug as well as 
the Majorana mass for the Ug. Above Yukawa couplings are automatically invariant under 
Lepton number as in Zee model. Lepton number as well as U{l)s symmetry is broken softly 
in Higgs sector through the following terms: 

fi^hth^ + %M2ht . (19) 

These soft symmetry breaking terms radiatively generate neutrino masses. The mass matrix 
generated at the 1-loop level has the following structure: 
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where 



ai 



( \ 

b a\ 

5 Al A2 

ai Al 5' 



ya2 A2 6' y 



C C 
fermi ; Al ~ —ff,r{ml - mi) 



(20) 



Vi Vi 

c 



(21) 



where we have neglected the electron mass. C and C in the above equations are given by 



C ~ ^f{mh) 



(22) 



/ and g are dimensionless functions of parameters in Higgs potential (denoted coUectively 
by m^) including the soft symmetry breaking terms displayed in eq. (p!9|) . 

Note that the 5, 5' are determined by the muon mass and rest by the tau as well as muon 
mass. Thus, eq.(^) has the perturbative structure required in our ansatz, eq.(H). One can 
calculate sphtting among neutrino masses generated by A^j, in eq. (PP]) . Treating 5,5' as 
perturbation and using the unperturbed eigenfunctions given in eq.(0), we get 

~ mi - {ócpsp + , 

~ -mi - [ÓCpSp + 5'CaSa) , 

m^z ~ ma + {òcpSp + ^'c^Sa) , 



my^ ~ -m2 + {5cpsp + 5'CaSa) ■ 



(23) 



This leads to the following sphttings: 



Aa ^ -Am2{5cpSp + 5'caSc) 
Ag _ mi 
Aa m2 



(24) 



The hierarchy in the solar and atmospheric scale is insensitive to strength of perturbation 
and is solely determined by the ratio of masses of the Dirac pairs. Eq.(pT|) leads to 

0-1 «2 fer 



Ai^ A2^ ff,r ■ 
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The required hierarchy in eq.(|T6D can be obtained if fer and /^^ are hierarchical. Eq.(|T 
now assumes the following form 



I ■ (25) 

'^1 A2' ferini (j^f^^TUr 



The difference on the LHS can thus be naturally small. It is possible to determine the basic 
parameters /e^^, /er, /^t and diretcly using eqs.f p^Sjpil) . Little algebra gives, 



C Vi Sa ' Pr mi SaCl 



_ m2Vi _ ^ c, 



fef, As mi 



mi C ' f^r AaSo, ' 

c„ . (26) 



Eqs.( p!5| , p4|) reveal that the required magnitudes of ai^2 niay be comparable to the per 



turbation 5, 5' . In this case, some of the perturbative results may change. It is therefore 
appropriate to perform exact diagonalization of matrix in eq.(^). We have done this nu- 
merically choosing values of parameters around the ones given in eq.(^). The specific values 
chosen are 

= -4.3 ■ 10-3; = 2.7 ■ 10-2 

/e^ = 2.1 ■ 10-^ /e, = 3.6 ■ 10-^; f^r = 8.6 ■ 10-6 • (27) 



We chose C ^ C ^ 0.01 for defìniteness. The above choice when substituted in eq.(pOD 
leads to the following values for the solar and atmospheric scales: 

A5 = 2.3 ■ 10-^ eV^ , 
A^ = 3.1 ■ 10-3 eV^ , 
Ai5^z5 = 0.3eV2 . (28) 



The mixing matrix is given by 



0.826 0.57 0.02 0.017 



-0.024 -0.011 0.71 0.71 



0.52 - .75 0.24 -0.2^ 



y -0.23 0.33 .64 -0.65 y 



(29) 



This can be seen to correspond to the following values for various parameters: 



sin^ = 0.99; sin^ 263 = 0.88 

sin^ 2eLSND = 2.7 ■ 10"3; sin^ 2e,hooz = 1-4 ■ lO'^ (30) 

si ^ 0.15 

The solar mixing angle is reduced signifìcantly compared to its maximal vaine in the absence 
of perturbation. This is welcome since strictly maximal mixing is not favored at least in the 
two generation analysis of the solar data [Q. Clearly, there would be ranges in the basic 
parameters of the model which would reproduces the correct mixing and masses. 

We note that the structure of the mass matrix in eq. (^0|) coincides with the one discussed 
in [|l^ but the underlying model presented here is much simpler. More importantly, phe- 



nomenological emphasis here is very different. We have shown that the basic structure in 
eq. (|20D displaying approximate Lg + — L,- — Lg symmetry provi des a concrete realization 
of the generalized 2 + 2 model which is fully consistent with ali the neutrino anomalies even 
after inclusion of SNO results. 



IV. SEESAW MODEL 



We now discuss how the ansatz of section (2) can be derived in seesaw type scheme. 
There are two ways of obtaining a light sterile neutrino in seesaw model. One is to assume 



that the mass matrix of the right handed neutrino is singular |T^. This can be done through 



some symmetry [^. The massless RH neutrino resulting from this singular matrix picks 
up a mass through its Dirac coupling with the active neutrino. The RH neutrino remains 
strictly massless and can provide the sterile state if its Dirac coupling is also forbidden by a 



symmetry. Example of this was recently presented in Alternative possibility is to add 
a sterile state to the conventional seesaw picture and impose symmetry which keeps it light. 
Examples of this possibihty were discussed in |]TD|,|T2[. Consider an active state ul, its RH 
partner upt and a sterile (left-handed) state with the following mass matrix in the basis 



( \ 





mcj 



(31) 
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This mass matrix leads to a massless, a light (~ — ^) and a heavy (~ M) neutrino. 
The sterile state mixes with the active state and influences the phenomenology. Cruciai 
points to note are the absence of Dirac mass term between and ul and the absence of 
the Majorana mass term for z/^. This can be achieved by means of some symmetry, e.g. R 
symmetry as in |TD[. This matrix was proposed as a model for solving the solar neutrino 



anomaly through — Ug mixing. Now we generalize the above idea to obtain the ansatz 
discussed in the last section. 

We consider the conventional seesaw picture with three left-handed and three right- 
handed neutrinos and add to it a sterile state Ug which would remain light. The lightness 
can be ensured by a structure which is generalization of eq.(PT|) to three generations. We 
demand separate conservation of Lg + — and the lepton number Lg corresponding to 
sterile state Ug. This would lead to Le + — Lr — Lg symmetric 4X4 matrix of the ansatz, 
eq.(^). We however need to break this symmetry softly in order to obtain realistic mass 
spectrum. Our model thus has two sets of mass terms: Cm which respect the symmetry and 
C'm which break it softly. The symmetric part is given by 

+ m^,^ Uf,L V^R + m^r VrL ^tR + H-C 

+ ^ {M,,vl^v,n + M,,vl^v,R + M^^vl^v^n + M^^v^^v^j,) . (32) 

Due to the combined effects of two symmetries, the Vg remains massless and decoupled 
from rest of the fermions in eq. (p2D . Its couplings with active neutrinos are induced entirely 
by the soft symmetry breaking sector CJ m chosen as follows: 

- = Pe^s^eR + Pt^i^s^tMR + Pr^gVrR + H.C. 

+ ]^M,eylRVeR . (33) 

Note that ali terms in the above equation connect only the sterile states and are soft in the 
technical sense- they would not lead to any divergences 0. Arbitrary choice of soft terms is 
consistent with this technical requirement. We have added ali possible singlet mass terms 
above except for the direct mass term v'gVg and some additional mass terms among the RH 
neutrinos. Addition of the latter terms do not make any qualitative change. The omission 
of v'gVg can be justified if soft terms are assumed to come through normalizable couplings of 



^ Models using such soft symmetry breaking were proposed in |22| 



11 



sterile state with standard model singlets r]i,r]2,ri carrying the (Lg + — Lr,Ls) charges 
(— 1, 1), (1, 1) and (—2, 0) respectively. Tliese fields cannot couple to u^Us but would lead to 
rest of the soft terms displayed in eq. (|5BD . In the following, we assume that Dirac masses 
Pe,ti,T are smaller than the scale in M^. This may be achieved through R symmetry as in 
T0| . Here we simply make this choice which is stable against radiative corrections. 

The effective 4X4 matrix emerging after seesaw mechanism can be written in the form 



of eq-d) with 



m 



M, 



/IT 



nie^nirr me^^iPr + PeZ + P^,Z^ / v) + rUeeiPeV + PfiZ) 
y^m^^rrirr m^^{pr + PeZ + p^^Z^ /y) + nif^eiPeV + PfiZ) j 



(34) 



yM^r 



{meeV + rrief^zf {rrieey + ruef^z) {m^eU + rn^^z] 

{rrieey + mef,z) {m^eV + m^^z) {m^eV + m^^zf 



\ 

ì 



(35) 



A' 







rrirrPpL "^PrPtM + 



rrirrPti 



Here, y = ^ and z = 

' ^ Alee Afe e 



Comparing eq.(|T^ with eq.(|3^ we fìnd that 

ai a2 m 



e/1 



Al A2 



< 1 . 



m 



provided the elements of the Dirac neutrino mass matrix obey the hierarchy 



e/1 



fili ■ 



(36) 



(37) 



(38) 



This hierarchy is a naturai assumption in many seesaw models. It would foUow, for example,if 
the Dirac mass matrix has Fritzch type structure. This structure also leads to 

ai 0-2 me^jn^e Pey + p^iZ 



Al A2 



Pr+PeZ+P^.z'^/y 



m. 



/i/i 



This hierarchy is dose to the one required on the phenomenological grounds, see 
eqs.(|^,|l3). Thus basic Le + — L^- — Lg symmetric ansatz is reproduced in this model 
under naturai assumptions. 

The symmetry breaking parameters Pe,fi,T are stili arbitrary. While many choices are 
possible, let us identify a particularly simple one. This corresponds to choosing z.p^j, much 
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less than the rest. The matrices A, A' of eqs. (|35|j36|) assume a simple form in the hmit 
mee,Pf_i, z tending to zero: 



A 



M, 



Ì1.T 







rnl^ ) 



A' 



M, 



/IT 







Vo pI ì 



(39) 



Treating non-zero elements in this matrix as perturbation we can evaluate the mass 
sphtting which turns out to be 

y 



As 
Aa 



-2mi 



y 



2 2 



UT 



( 2 2, 2 2 \ 



jjLe 



m2c}pl + ci...^^ 



^2™ 2 



(40) 



This can be consistent with the phenomenological requirement for proper choice of pa- 
rameters. 

The above exercise though illustrative in nature, shows that it is indeed possible to 
integrate a singlet neutrino in seesaw picture in a way that leads to generalize 2 + 2 model 
and particular ansatz discussed in section(2). Gross features of the ansatz, particularly 
the required Le + — Lr — Lg symmetric matrix m follows under the standard seesaw 
assumptions. The full matrix with a broken L^. + — ~ Lg symmetry is complex but 
proper symmetry breaking can be achieved due to large number of parameters in the model. 



V. SUMMARY 



The three sets of experimental results namely, the observed solar and atmospheric neu- 
trino fluxes as well as probable oscillations seen at LSND calls for a consistent theoretical 
explanation. 2 + 2 schemes with four neutrinos were considered an attractive mechanism to 
explain ali anomalies. In these schemes, the z/g or i/^ was exclusively assumed to convert to a 
sterile state. This scenario is not supported by the results of the solar neutrino experiment 
at SNO. Simple generalization which can stili explain ali the experimental results assumes 
that sterile neutrino flux is simultaneously but partially present in the solar as well as atmo- 
spheric neutrino fluxes. We have shown that this possibility follows naturally when neutrino 
mass matrix is + — Lr — Lg symmetric. Assumption of this symmetry in fact allows 
US to completely reconstruct neutrino mass matrix from the experimental results. We have 
undertaken this exercise in this paper. We have also shown that resulting structure follows 
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naturally in two of the conventional schemes based on Zee model and the seesaw model for 
neutrino masses. 

Note:Before posting this paper, we noticed a preprint ( |hep-ph /Oli 0243] ) by K.S. Babu and 
R. N. Mohapatra who also advocate the use of Le + — — Lg symmetry as an explana- 
tion of the neutrino anomalies. Model presented in this paper is different from the models 
considered here. 
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